EPJ manuscript No. 

(will be inserted by the editor) 



Social dynamics with peer support on heterogeneous networks 



The "mafia model" 



Marta Balbas Gambra and Erwin Frey 



o 

< 
(N 



43 

6 
o 

O 

OH 



> 

in 
oo 
in 

o 



X 



Arnold Sommerfeld Center for Theoretical Physics and CeNS, Department of Physics, 
Ludwig-Maximilians-Universitat Miinchen, Theresienstrasse 37, D-80333 Miinchen, Germany 

Received: date / Revised version: date 

Abstract Human behavior often exhibit a scheme in which individuals adopt indifferent, neutral, or radical 
positions on a given topic. The mechanisms leading to community formation are strongly related with 
social pressure and the topology of the contact network. Here, we discuss an approach to model social 
behavior which accounts for the protection by alike peers proportional to their relative abundance in the 
closest neighborhood. We explore the ensuing non-linear dynamics emphasizing the role of the specific 
structure of the social network, modeled by scale-free graphs. We find that both coexistence of opinions 
and consensus on the default position are possible stationary states of the model. In particular, we show 
how these states critically depend on the heterogeneity of the social network and the specific distribution 
of external control elements. 
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1 Introduction 

Opinion and group formation in human societies have 
been investigated in the last years as complex phenom- 
ena, well described with the methods of non-linear dy- 
namics and statistical physics (for a review see [I]). It has 
become clear that the evolutionary dynamics of societies, 
based on social pressure and imitation, must be under- 
stood together with the intricate network of contacts es- 
tablished among single agents. Individuals are represented 
by the nodes of a graph and interact with their neighbors 
under a given set of rules to form their own opinion. Inter- 
actions may take place in several ways, and many models 
have dealt with various mechanisms in which individuals 
might meet as well as with different rules to update their 
opinions. 

In the simplest voter model [21GIH individuals on the 
nodes of a regular lattice may hold two opinions encoded 
by a binary spin variable a = ±1. Two randomly selected 
individuals interact with each other (one-to-one interac- 
tion) in which the second simply adopts the opinion of 
the first. In one dimension, the voter model is identical 
to the kinetic Ising model with zero-temperature Glauber 
dynamics (SJIHIEI- The latter is exactly solvable for reg- 
ular lattices in any dimension [8,9J. This simple dynam- 
ics leads to consensus in finite systems, which is the only 
stable solution: all spins are aligned, while both states are 
equally likely to be the consensus opinion — provided they 
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are equally abundant in the initial state — i.e. ensemble 
magnetization is conserved [5] . 

On heterogeneous networks, special updating mecha- 
nisms like the link-update [TUlllll[T2"] approach are needed 
to ensure conservation of the ensemble magnetization. For 
the reverse-voter [71lll|. where a node is randomly selected 
and its opinion copied to a randomly selected neighbor, it 
was found that the time to consensus differs from the dir- 
ect voter model. A common feature of the voter dynamics 
on various complex networks, such as in small- world |13| , 
scale- free |10l[Tn[T2] . and random [T5] graphs, is the ab- 
sence of consensus in the thermodynamic limit, which has 
also been supported by analytical investigations on uncor- 
related graphs [TS]. However, consensus is systematically 
reached on finite graphs [121113] . even though metastable 
states of coexisting opinions are observed for long periods 
before the system evolves towards consensus for a time 
which scales with the system's size. 

In the voter model, the inclusion of more than two 
opinions [161117] , non-confident vacillating voters [TS] , zea- 
lots with a fanatic position [T91I20], or a threshold num- 
ber of successful encounters [2T| allows the coexistence of 
opinions. Encounters among many agents have been con- 
sidered in the majority model [221123] . in which a group 
of individuals is randomly selected and all of them ad- 
opt the opinion of the majority (all-to-all interactions), 
or in the rumor spreading model [241125] . in which more 
than one group-encounters take place simultaneously. In 
both of these models consensus on one or the other opi- 
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nion depending on the initial conditions is achieved for 
various structures. Nevertheless, if noise is introduced in 
this group dynamics, as in the majority-minority |26| or 
the majority-vote [5] models, the system evolves either 
towards a consensus or towards a state in which both opi- 
nions are equally represented (zero magnetization state in 
the Ising model analogy). The same behavior is observed 
in the Sznajd model |27U28| where two individuals sharing 
the same opinion impose it on their neighbors. 

Another approach for describing social interactions that 
accounts for the influence of the whole neighborhood on 
a single individual {one-to-all interactions) has been dis- 
cussed in social impact theory [2DII301I3T]. An agent changes 
her opinion if the pressure in favour of the opinion change 
overcomes the support to keep the current position. This 
model yields stable coexistence of opinions, unless the 
presence of external fields is taken into account, for which 
metastable states are observed in which domains with the 
minority opinion successively shrink [3 1 J . 

The Abrams-Strogatz (AS) model describes the com- 
petition between two languages in a population, with non- 
linear transition rates proportional to x a , where x is the 
population fraction of one particular language, and a the 
volatility characterizing the tendency to change state [32] , 
For high volatilities, a < 1, coexistence of both languages 
has been found, whereas for a > 1 one of the two lan- 
guages becomes dominant [35] . Studying the AS model 
on networks, it has been shown that a decrease in network 
connectivity leads to a reduced parameter regime with lan- 
guage coexistence [33 . In the same study, the AS model 
was also compared with a variation including a bilingual 
state, which has been observed to hinder coexistence. 

Here, we propose a model which aims to capture the 
mechanisms of community formation and how these de- 
pend on the specific structure of the contact network. It 
describes the dynamics of a society in which individuals 
may either be neutral (hold a default opinion) or belong 
to a minority with radical positions on a current topic. 
Examples of human behavior with such a dynamics are 
the membership and involvement in political parties or 
organizations, belonging to religious communities, the in- 
terest in leisure activities as online games, the behavior of 
consumers with new technologies, or the typical cycle of 
addiction to chemicals. Neutral individuals may at some 
point in time make acquaintance with a new product or 
idea and subsequently become customers or adopt it as 
their default position. Hereafter, they can either like it 
and consolidate as their default state or not and return to 
neutral behavior. Additionally, if agents in a default state 
are in contact with fanatics (radical agents) they can also 
turn into radicals. Radical agents, in turn, may get bored 
or disappointed after a while reducing their level of com- 
mitment and turn back to a default position or even com- 
pletely leave the community and become neutral agents 
again. 

There are two characteristic features of our model, 
defined mathematically in the next chapter, which are es- 
sential for its dynamics. First, we account for the fact that 
the population density may vary over the contact net- 



work, i.e. in addition to individuals with different states 
(opinions) we also consider empty sites. Second, the prob- 
ability of a given individual to change its state is taken 
to be proportional to both the number of like and alike 
peers in its neighborhood. This product form of the trans- 
ition rates is the main difference to previously studied non- 
linear voter-like models [?JM] and the AS model [32]. This 
form implies, that transition rates may change from a lin- 
ear to a non-linear form upon changing the composition 
of the population. 

The outlined scheme bears similarities with the form- 
ation of mafias, in which normal citizens are the analog to 
default opinions, while mafias represent the radical minor- 
ity. Neutral individuals are those who do not take active 
part in the society and that are represented by empty 
sites. Allegorically these sites may enter the community 
as citizens by a birth event and join the mafia later due 
to the social pressure of their vicinity. In this paper, for 
illustration we will use the metaphoric language of mafias. 

We will show that the mafia model yields a rich phase 
behavior which contains both coexistence of opinions and 
consensus on the default position in well-mixed societies. 
In scale-free networks, coexistence appears in a larger re- 
gime of the parameter space due to the effect of local in- 
teractions: the stability of the absorbing state (extinction) 
is lost. We will discuss the role played by the nature of 
interactions — one-to-all, i.e. the whole neighborhood de- 
termines the state of single agents, as compared to one- 
to-one interactions like in the voter model — for the loss 
of stability of the absorbing state. Furthermore, we will 
show that the inclusion of external fields in the model, 
representing the presence of elements controlling the in- 
teractions (which we refer to as police individuals), gives 
rise to micro-phase separation where regions with different 
dynamics are observed to coexist in the system. 

This paper is organized as follows: We first describe 
the dynamics of the mafia model and then we analytically 
discuss the phase diagram expected for well-mixed popu- 
lations. Then, we analyze the behavior on structured soci- 
eties showing different degrees of heterogeneity, which we 
investigate through extensive stochastic simulations. We 
finally address the role of external control elements in the 
system and explore how specific distributions of those into 
highly or sparsely connected nodes modifies the system's 
behavior. An extended mean-field theory accounting for 
the local structure of networks is given in the Appendix. 

2 The mafia model 

We model human interactions in social networks by graphs 
whose nodes allocate interacting agents and whose edges 
describe the relations between pairs of individuals (Fig.[T]). 
The nodes of a network can either be empty (<P) or occu- 
pied by an individual. In the simplest case agents belong to 
one of two existing groups identified with what we refer 
to as two different strategies, namely they either belong 
to the mafia (M) or to the group of lawful citizens (C). 
The network's size is assumed to be fixed and, therefore, 
the sum of the number of citizens, mafiosi, and empty 
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Figure 1. [color online] Individuals (red mafiosi, green cit- 
izens, white empty sites) are distributed in the nodes of a 
scale-free network. The composition of the neighborhood of 
the agent who makes a decision, marked with a question mark, 
reads (c, m, 4>) = (3/5, 1/5, 1/5). 

places (or their corresponding fractions) remains constant: 
C + M + <P = N. 

We employ an agent-based model in which the society 
evolves according to the following stochastic dynamics: (i) 
at empty nodes citizens are born with a site-independent 
fixed rate b, (ii) at occupied nodes individuals die with a 
constant rate d, and (iii) strategy changes take place in 
which citizens join the mafia at rate w c -y m {=: w cm ) while 
mafiosi leave it at a rate w m -> c (=: w mc ). Since individu- 
als who reconsider their status, namely whether leaving 
or joining the mafia, are influenced by their surround- 
ings, these rates are determined by the composition of 
the neighborhood of a given individual. 
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As the neighborhood of a given agent we define the set 
of all its nearest neighbors, i.e. the set of all other agents 
connected directly to it through an edge of the network. 
In Fig. [TJ for instance, the neighborhood of the citizen 
with the question mark contains three citizens (green), 
one mafioso (red), and one empty place (white). The com- 
position of this neighborhood is characterized by the tuple 
(c, m, </>), where c, to and <p denote the fraction of mafiosi, 
citizens and empty sites in the neighborhood of a given 
individual at a given instant of time. 

Single agents feel both the pressure from neighbors 
holding a different strategy and the support of neighbors 
who are members of their peer group, motivating the fol- 
lowing choice for the transition rates: 

Wc-ym = w cm = m s m (1 - c), (1) 
Wm^-c = w mc = cs c (1 — to). (2) 

In this way, the pressure exerted on a given agent by 
individuals belonging to the adversary group is propor- 
tional to their fraction among the nearest neighbors, i.e. 



to and c, respectively. The strength of this pressure, s m 
or s c , which we also call persuasiveness, is taken to be 
the same and constant for all individuals. At the same 
time, the supportiveness of the alike individuals attenu- 
ates this pressure through a multiplicative factor, (1 — c) 
and (1 — m), respectively. This product form of the trans- 
ition rates implies that individuals do not reach a stable 
state but keep on changing their strategies unless the full 
neighborhood holds the same opinion and the agent is 
fully surrounded by alike agents. As a consequence the 
relative strength of the linear and non-linear term in the 
transition rates change with species frequencies, cf. uj cm ~ 
to(1 — c) = m(m + 4>): If the abundance of the adversary 
species is much smaller than the fraction of empty places 
(to <C 4>)i the leading term to change an indiviuals' state 
is linear, whereas for species abundances significantly lar- 
ger than the fraction of empty places the leading term is 
non-linear and the dynamics changes accordingly. This is 
distinct from previously studied non-linear opinion mo- 
dels like the Abrams-Strogatz model where the type of 
non-linearity remains fixed independent of the composi- 
tion of the populations [32, 33J . In our notation the rates 
of the AS model read u> cm = s m m a and uj cm = s c c a with 
volatility a, and prestiges s m / c of "languages" m and c. 

Moreover, since the tuple (c, m, (f>) denotes the over- 
all composition of an individual's neighborhood, a given 
agent interacts with its neighbors in a one-to-all scheme. 
As it will turn out, the ensuing dynamics and stationary 
states on networks are genuinely different from those ob- 
tained for a one-to-one scheme where a given individual 
interacts pairwise with a single randomly chosen neighbor- 
ing site. The main difference is that in one-to-all schemes 
one has to account for a discrete set of possible neigh- 
borhood compositions. This feature invalidates standard 
mean-field theories even in models where one adds strong 
mixing. As shown in Appendix A, an extended mean-field 
theory accounting for these discrete set of possible neigh- 
borhood compositions explains why a one-to-all interac- 
tion in structured societies induces changes in the phase 
diagram. 

In this paper, we will focus on the asymmetric case in 
which mafiosi are more persuasive (much stronger) than 
citizens, s m 3> s c ps 0. For this case, the transition rate 
w mc vanishes and individuals can leave the mafia only 
rather indirectly, namely via a death-birth process. The 
more general case yields less interesting dynamics; a full 
discussion may be found in |35j . 

Decision processes may, in general, also be influenced 
by many external factors, such as mass media, independ- 
ent of the interacting actors. More particularly, societies 
may provide means to regulate their conflicts and protect 
citizens against the damage of mafias. One would like to 
account thus for the role of control elements such as po- 
lice in actual societies. Control elements do not particip- 
ate directly in the social dynamics by interacting with the 
agents and evolving in time, but specifically control the 
relations between pairs of individuals. They are located 
at the edges connecting nodes between two individuals of 
a social network. The fraction of total edges allocating 
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control elements is p. The model includes their two-fold 
catalytic role activating the transition m — > c and inhib- 
iting the transition c —> m: 

w cm =ms m (l-c)(l-p), (3) 

W m .c = (CS C + S p p)(l - 17l). (4) 

The control elements protect citizens from the mafia's 
pressure, attenuating their strength proportional to the 
factor (1 — p). Simultaneously they persuade individuals 
to leave the mafia proportional to their fraction p and 
strength s p . The persuasiveness adds to that of the cit- 
izens, both being still attenuated by the presence of mafi- 
osi as before. 



3 Mean-field approximation 

For well-mixed populations in which every individual in- 
teracts with all other agents in the system, the problem we 
pursue to solve is topologically equivalent to having the 
agents allocated on the nodes of a complete graph. Then 
the local field which every agent experiences is the same 
as the global field, given by the average relative abund- 
ance in the whole population. This simplification leads 
to a mean-field approach for well-mixed societies which 
allows an analytical characterization of the problem. We 
note that this type of mean-field limit does not correspond 
to the limit of fast diffusion on networks if the interaction 
is one-to-all as is the case in our mafia model. Then, the 
dynamics of the network-structured society is actually bet- 
ter described within an extended mean-field theory which 
accounts for the fact that for each given node there is a fi- 
nite set of neighborhood compositions, cf. Appendix A. In 
the following we discuss the standard mean-field approach 
in order to highlight the novel effects introduced by the 
network structure and the one-to-all interaction scheme. 

We non-dimensionalize by choosing the inverse death 
rate, d , as our basic time scale and define dimensionless 
time and all other parameters accordingly: r = td, (3 = 
b/d, o~i = Si/d, Uij = Wij/d, with i,j £ {c, m}. Then, the 
mean-field equations for the time evolution of the different 
population fractions read: 

c= (f>/3 + mu) mc - cu) cm - c, (5) 
m = -mu mc + cw cm - m, (6) 
= -00 +(1-0, (7) 

where a dot signifies a time derivative with respect to r. 
From Eq. ([7]) we immediately see that the stationary frac- 
tion of empty places is independent of the social dynamics 
encoded in the transition rates u cm and u) mC) but only de- 
pends on the dimensionless birth rate /3: 



In other words, the birth rates determines how densely 
populated the system is. Inserting this result together with 
the constraint c + m + (j) = 1 into Eq. © yields an implicit 



equation, whose solutions are the fixed points or station- 
ary states of the population dynamics: 

muj mc (c,m,p) - ^Y~4~g - m ) w cm (c, m,p) + m = 0. 

In the following, as noted above, we focus on the asym- 
metric case where mafiosi are much stronger than citizens, 
er rn 3> o~ c ~ 0, and use the simplified notation tr := a m . 
Then, in the absence of control elements, the mean-field 
equations reduce to: 

c = 4> j3 — am (I — c) c — c, (8) 
m = a m (I — c) c — m, (9) 

whose fixed points are: 

m° = 0, m ± = I(_^-l ±v /i— 47^). (io) 

The ensuing phase diagram (Fig. [5]) shows three differ- 
ent regimes: coexistence, bistability, and mafia extinction. 
For a given birth rate /3, mafiosi unavoidably get extinct 
below the threshold strength cr* n , while for strengths lar- 
ger than <j* c coexistence is the only stable solution. In the 
intermediate regime, a* n < a < a% c , the system is bistable 
and the stationary solution critically depends on the ini- 
tial conditions. It turns out that the minimal strength of 
mafiosi to ensure that they do not get extinct increases 
with the birth rate of citizens, although very weak mafiosi 
die out independently of the birth rate. 
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/? 

Figure 2. Phase diagram as a function of the birth rate j3 and 
persuasiveness a. Three regions with different regimes (coex- 
istence, bistability, extinction) are shown, separated by a tran- 
scritical (solid) and a saddle-node (dashed) bifurcation. 

For a particular value of the birth rate (3^0 (Fig. [3]) , 
the absorbing state m° becomes unstable via a transcrit- 
ical bifurcation for a persuasiveness larger than the threshold 
value <7* c = (1 + /3) 2 //3, while the coexistence state m + is 
real only above the saddle-node bifurcation at a* n = 4 
(for (3 > 1). In the regime between these two values the 
system is bistable. The third solution m~ is unstable for 
all meaningful values. The stationary fraction of mafiosi 
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in the coexistence state increases with their strength up 
to the asymptotic value m + =p = /3/(l + /3)in the limit 
of infinite strength, a —> oo, where p is the fraction of 
populated nodes. 
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Figure 3. Bifurcation diagram for a fixed birth rate /3 — 10 as 
a function of the strength parameter a. The solid and dashed 
lines represent stable and unstable solutions, respectively. The 
two solid circles stand for the saddle-node bifurcation, at which 
the solutions m become complex, and the transcritical bifurc- 
ation, above which the absorbing state m° becomes unstable, 
respectively. In the dark grey and white areas the system is 
monostable — extinction and coexistence are the stable states 
respectively — and bistable in the light grey region. 



4 Structured social networks 

Individuals in actual societies do not interact with the 
whole population, but rather only with agents in their im- 
mediate vicinity. A given agent who updates her strategy 
is influenced by her nearest neighbors, so that the frac- 
tion of citizens, mafiosi, and control elements entering the 
transition rates ([!} and ^ are now given by the instant- 
aneous composition of her neighborhood. The population 
frequencies (c, m, cf>) are restricted to a, finite set of possible 
combinations determined by the size of the neighborhood. 
They are thus no longer continuous as in the mean-field 
case discussed in the previous section, but take a discrete 
set of values determined by the composition of the neigh- 
borhood of a given node. 

We investigate the behavior of the mafia model on 
scale- free networks (SFN) whose degree distribution scales 
as a power law p(k) oc fc -7 . The parameter 7 is related to 
the heterogeneity of the network and increases for more 
homogeneous structures. Scale- free networks [36 , 37J ex- 
hibit very interesting features suitable for the modelling 
of social systems, such as large heterogeneity — with a non- 
negligible fraction of highly connected nodes — and relat- 
ive small average path-lengths, i.e. the average over the 
shortest distance between all pairs of nodes. The latter 
scales logarithmically with the network's size, as it has 
been observed in human communities |36j . 



The scale-free networks used in this work have been 
generated following the uncorrelated model [38] , for which 
the resulting average degree (k) is a function of 7. This 
method imposes a maximal cutoff for the degree, auto- 
matically yielding the average degree for a given network 
size; see [35] for further details. In this work, we have used 
networks with increasing 7 ranging from 2 to 4, and cor- 
responding decreasing average connectivities (k) varying 
from 7.6 to 2.4. 

We have performed stochastic simulations employing 
an agent-based model following the dynamics described 
above. Individuals belonging to both populations were 
taken as randomly distributed with specific initial frac- 
tions mo and cq. We used random sequential updating ac- 
cording to the dynamical rules defined previously. The res- 
ults discussed in this paper were averaged over 1000 runs. 
Starting from a given initial state, the system evolved un- 
til it reached a quasi-stationary state; since the system has 
an absorbing state, m° = 0, it will be unavoidably reached 
if one waits long enough. The population remained in this 
state for a long enough time window where the relevant 
observables were measured. We have carried out a finite 
size analysis and assessed that the systems reaches a quasi- 
stationary state at a time r' proportional to the system's 
size, t' = XN, with the proportionality factor A ranging 
from 0.01 to 0.015 (not shown). 

We characterize the stationary state of the system in 
terms of the extinction probability p CK t, which is the prob- 
ability that mafiosi got extinct at some time r' within 
the quasi-stationary time window. The ensuing phase dia- 
grams, as shown in Fig. 0J are qualitatively different for 
well-mixed graphs and scale-free networks. In particular, 
the topology of the underlying network strongly matters 
in the parameter regime where mean-field theory predicts 
bistability, i.e. a dependence of the stationary state on 
the initial state of the population. We have specifically 
analyzed the limiting case in which the initial population 
of mafiosi vanishes (mo — > 0). In this limit, the stand- 
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Figure 4. [color online] Extinction probability — as encoded in 
the side bar — for initial conditions with a very small fraction 
of mafiosi, mo = 5 ■ 1CF 4 , for (a) a complete graph of size 
N — 6000 and (b) a scale-free network (SFN) with 7 = 3 and 
N = 8000 sites ((k) = 3.16). The solid black line represents 
the deterministic mean-field prediction separating a region of 
coexistence above it from a regime in which mafiosi get extinct 
below it. 
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ard mean-field theory asserts that mafiosi go extinct for 
all of the bistable area. This is indeed the case for well- 
mixed societies (complete graphs), where the results of 
our stochastic simulations closely resemble the prediction 
of the (deterministic) mean-field theory. There are some 
minor differences in the actual position of the phase bound- 
ary which we attribute to stochastic effects; see Fig. Sk- 
in contrast, coexistence is systematically found in the bi- 
stable regime (of mean-field theory) for all studied net- 
works (scale-free networks with 7 ranging from 2 to 4)0- 
An exemplary phase diagram for a network with 7 = 3 
and a small initial mafiosi fraction is shown in Fig. \^p. 

We attribute this interesting feature of what one could 
call an "evolutionary stability" of mafiosi to the locality of 
interactions on networks or regular graphs. In structured 
societies, the one-to-all character of interactions becomes 
relevant, as averaged frequencies suitably reproduce inter- 
actions between two individuals (one-to-one) but fail to 
describe interactions with a whole vicinity. An extended 
mean-field analysis, accounting for the discrete set of pos- 
sible neighborhood compositions is thoroughly discussed 
in Appendix A. This extended mean-field theory gives a 
phase diagram where the bistable regime is drastically re- 
duced in favour of the coexistence regime, in accordance 
with the above numerical results. 

What is the micro-dynamics giving rise to these res- 
ults? A vanishing small fraction of mafiosi gives rise to neg- 
ligible transition rates iu cm within a standard mean-field 
approach and consequently the mafia dies out in a short 
period of time. However, if interactions are local, citizens' 
vicinities containing any mafioso yield considerable trans- 
ition rates for updating citizens to become mafiosi. Under 
which circumstances can a few isolated mafiosi invade a 
citizen population? The following discussion gives a set of 
heuristic arguments which are not meant to be quantitat- 
ive but illustrate the basic mechanism behind the observed 
evolutionary stability. 

Starting from some isolated mafiosi, the dynamics may 
be separated into three phases as illustrated in Fig. [SJ 
invasion-expansion-saturation. First, the isolated mafiosi 
must invade some neighboring cell before dying, i.e. the 
probability for any citizen in the mafioso's neighborhood 
to change her strategy must be larger than her death prob- 
ability. The probability for a mafioso at site i to invade 
any neighboring cell is the sum over the probabilities for 
all her ki neighbors (assumed to be citizens) to change 

strategy, i.e. p inv = J2jLoPi->-m = J2f=Q (jm ( 1 ~ c ) At i 
which yields kia(l/kj) 2 At considering that every citizen 
has a single mafioso in her vicinity. If the structure is 
homogeneous enough, so that one may consider similar 
connectivities, i.e. ki ~ A;., then the invasion probability 
becomes a /ki At. The invasion probability is larger than 
the death probability for mafia's strength larger than the 
degree cr t i nv > fc,. This heuristic argument nicely explains 
the results obtained from simulations which show that the 

1 We have also performed simulations on square grids which 
show similar behavior as scale- free networks |35] . 

2 Note that this approximation is formally valid only for SFN 
with large 7 or random graphs. 
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Figure 5. [color online] Invasion-expansion-saturation process. 
Invasion: A few isolated mafiosi invade their nearest neighbors 
instead of dying, which leads to the formation of very small 
clusters of mafiosi. Expansion: these small clusters expand as 
more and more citizens successively join the mafia. Saturation: 
the society reaches a stationary state with a dynamic balance 
between birth and death processes. 

extinction probability decreases with the network's av- 
erage degree (k) (not shown). Say now the mafioso has 
not died out but managed to invade some neighboring 
cells and thereby created small mafia clusters. In order for 
these clusters to survive they have to expand as groups to 
achieve a stable coexistence state in the population. Mafia 
clusters will grow if the fraction of citizens becoming mafi- 
osi at their interfaces is larger than the fraction of those 
citizens gained via birth process. Assuming similar frac- 
tions of citizens and mafiosi in the vicinity of an agent loc- 
ated at the interface, m ~ c~ 1/2, and large enough birth 
rateifl/3, the fraction of born citizens is given by li n t4>P and 
those becoming mafiosi /i n tw cm ~ l ln tcrm(l — c) ~ hnt& /4, 
where k at is the length of the cluster interface. Mafia 
clusters will thus grow if their strength becomes larger 
than some threshold value, a > Ctcxp = 4. Finally, mafia 
clusters will expand to the point in which there are no 
citizens left at which expense they can grow. This satura- 
tion takes places for a residual fraction of citizens, which 
equals the fraction of born citizens at empty places. A 
dynamic stationary state has thus been reached at which 

3 This assumption implies that one can neglect empty sites 
and approximate the fraction of citizens and mafiosi at both 
sides of the cluster boundary by 1. The approximation becomes 
better with an increasing size of the vicinities. For small neigh- 
borhoods, taking empty sites into account would make it easier 
for the mafiosi to expand their territory as citizens have less 
peer protection. Therefore, accounting for empty sites would 
slightly lower the threshold strength such that our approxim- 
ation can be regarded as an upper bound. 
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both populations coexist due to the asymmetric birth of 
citizens — otherwise citizens would get extinct. These heur- 
istic arguments for the growth process in structures, sup- 
ported by the results of the stochastic simulations, nicely 
evidences that above some threshold value for the mafia's 
strength, the mafia can always invade a citizens popula- 
tion. This threshold is given by the larger of the following 
two values: the mafia strength required to invade neigh- 
boring sites, Otinv) or for small clusters to expand, Otoxp, 
i.e. a > max(0-ti nv ,crtexp) = max( (fc) , 4) . 

In the above line of reasoning, we have assumed that 
very node has the same connectivity, which is true only 
for random networks. However, since mafia's expansion 
is a process which depends on the local dynamics, we 
suppose that our heuristic approach should remain valid 
for scale-free networks. Moreover, the existence of regions 
with a connectivity significantly smaller than the average 
in scale-free networks suggests that coexistence is more 
likely to be found there than for the corresponding ran- 
dom graph with the same average degree. 

Another interesting result is that in the limit of com- 
plete random graphs, (k) — > N, the above heuristic ap- 
proach gives a threshold strength a > 1 + /? for mafiosi 
to survive, which agrees well with the standard mean-field 
analysis — see Fig. H^i. 



5 The role of control elements 

Societies can regulate themselves by introducing control 
elements which fight the expansion of mafias. Here, we 
consider control elements (police) which are randomly at- 
tached to the edges of a given graph, such that a frac- 
tion p of the total number of edges — and thus of pairs 
of individual relations — are controlled. With an increas- 
ing number of police in the asymmetric model the mafiosi 
also need to enhance their strength/persuasiveness in or- 
der to survive. For simplicity, we consider parameters such 
that a m = f3 = a corresponding to the diagonal in Fig. [2] 
where only two regimes are observed in the mean-field 
limit, namely extinction and bistability. 

A standard mean-field analysis - equivalent to a well- 
mixed population or alternatively a complete graph - in- 
cluding control elements with strength a p = a shows the 
same type of regimes as for the model without police: For 
small police fractions the system is bistable and coexist- 
ence or extinction is obtained depending on the initial 
composition of the population. Above a certain police frac- 
tion, the mean-field analysis predicts extinction of mafiosi 
as the only stable fixed point; see Fig. [5^ for a phase dia- 
gram with initial condition mo = 0.5. The mafia's strength 
required to achieve coexistence is increased with the pres- 
ence of control elements. 

On scale-free networks, the local character of interac- 
tions again leads to an enlarged coexistence region as com- 
pared to the well-mixed case, cf. Fig. [SJd. However, there 
is a significant difference to networks without control ele- 
ments. In the present case, the coexistence regime is en- 
larged beyond the bistable regime of mean-field theory and 
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(a) complete graph (b) SFN 7 = 4 



Figure 6. [color online] Extinction probability for a policed 
system on (a) a complete graph (well-mixed population) with 
N — 6000 and on (b) the nodes of a scale- free network with 7 = 
4 ((k) = 2.45) and N = 8000. The solid line in (a) represents 
the separatrix predicted by the deterministic mean-field theory 
for the given initial conditions mo = 1/2. 



now also covers an area in parameter space where the sys- 
tem was monostable for well-mixed populations (i.e. where 
extinction was the only stable state); cf. Fig. [6^1, and b. The 
reason for this anomalous behavior is that the control ele- 
ments introduce a kind of quenched disorder, which from 
thermodynamic systems is known to strongly affect phase 
behavior. Indeed, we have checked by numerical simula- 
tions (data not shown) that coexistence is found only in a 
parameter regime where mean-field theory predicts bista- 
bility, if the control elements are allowed to diffuse on the 
network (annealed disorder). Analyzing the population 
dynamics on a local scale, we observe that the presence 
of control elements leads to a kind of micro-phase separ- 
ation in the society: while citizens preferentially populate 
policed regions, mafiosi aggregate in unpoliced areas. The 
unpoliced areas display the same dynamics as the fully 
asymmetric model, in which species coexist for mafias' 
strength larger than some threshold, a > at- Fig. [7] shows 
a snapshot of a typical population structure on a square 
lattice, where both regimes are clearly recognizable. 




Figure 7. [color online] Snapshot of the population structure 
on a lattice with a police fraction p = 0.3 and size N = 4900. 
Red cells are occupied with mafiosi, green with citizens, and 
white cells correspond to empty nodes. The black segments 
represent policed connections between adjacent cells. Citizens 
and mafiosi cluster preferentially around policed and unpoliced 
areas, respectively. 
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Micro-phase separation results in an increase of the 
respective minority's population size with increasing net- 
work homogeneity 7, or, equivalently, decreasing average 
degree (k), cf. Fig. [5] This can be understood as follows: 
Increasing network homogeneity implies smaller neighbor- 
hood sizes and larger average paths between pairs of nodes, 
and thereby increasingly hinders mixing of mahosi and 
citizens. As a consequence, phase separation into isolated 
mafia and citizen clusters becomes more pronounced. Then, 
the respective minority species can find niches where they 
can grow to population sizes larger than in a correspond- 
ing well-mixed environment. 



a 



p — 0.1 mafiosi 
p = 0.1 citizens 




p = 0.2 mafiosi 
0.2 citizens 
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Figure 8. [color online] Stationary fractions of mafiosi and cit- 
izens as a function of the network homogeneity 7 for two differ- 
ent values of the police fraction p. The values of p were chosen 
such that for 7 = 2 the stationary state corresponds either 
to coexistence (p = 0.1, red), or mafia extinction (p = 0.2, 
blue). Hence, for p = 0.1 citizens are the minority species, 
while for p — 0.2 mafiosi are the minority. Solid and open 
symbols represent mafia and citizens fractions, respectively. In 
both cases the respective minority species increases in popu- 
lation size with increasing network homogeneity 7, or, equi- 
valently, decreasing average degree (k). The average degrees 
with increasing 7 ticks read (k) — 7.63,4.31,3.16,2.69,2.45. 
Data are obtained for the following parameters: N = 10000, 
m = c = 1/2, a = f3 = 10. 



A finite size analysis of the extinction probability with 
the increasing heterogeneity of graphs, for a fixed police 
fraction (p = 0.3), is shown in Fig. [9] The increasing 
slope of the extinction probability with N suggests that as 
N — > 00 there might be a phase transition separating an 
extinction from a coexistence regime. To make this point 
conclusive, however, would require to simulate even larger 
system sizes. Our main point is that for sufficiently homo- 
geneous networks (large 7) the phase separation between 
policed and unpoliced areas allows survival of mafiosi in 
the unpoliced areas for police fractions which otherwise 
would suffice for their extinction in well-mixed and het- 
erogeneous populations. 
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Figure 9. Extinction probability of mafiosi for a fixed po- 
lice fraction p = 0.3 as a function of network homogen- 
eity 7. With increasing overall population size N the trans- 
ition between certain extinction and certain survival of mafi- 
osi becomes very pronounced indicating that there might be 
a phase transition as a function of 7. Data were obtained 
for the following set of parameters: a — ft = 10 and mo = 
co = 1/2. The average degrees with increasing 7 ticks read 
(k) = 7.63,4.31,3.16,2.69,2.45. 



Targeted police distribution 

In the same way that vaccination during an epidemic is 
more effective if targeted to highly connected nodes or 
hubs |39U40| . one might suppose that the distribution of 
control elements also matters in eradication of mafias. To 
explore this question we investigate the effect of two spe- 
cific distributions q a {eij) of control elements: 



didj 



Eij d i d j ' 
(didj)- 1 



" A " J) V,7,/,/ ; : : 



(11) 
(12) 



These probability distributions for a control element to 
be allocated at an edge connecting nodes i and j are 
taken as either directly (qi) or inversely ((72) proportional 
to the product of the degrees of the nodes connected by 
the edge, didj. Therefore, the first distribution favours po- 
lice attachment close to highly connected nodes, whereas 
the second promotes control elements to surround sparsely 
connected nodes. 

The network topology should be the key element for 
understanding the relevance of a specific police distribu- 
tion. Following the analogy with the problem of vaccin- 
ation, one would expect that surrounding network hubs 
with control elements should hinder the proliferation of 
mafias, since individuals in hubs are expected to be more 
influential. Although highly connected individuals influ- 
ence the decision process of a larger number of agents, 
their influence on the targeted individual is as import- 
ant as that of sparsely connected agents. Therefore, con- 
trolling the population in hubs is not decisive for the sys- 
tem's dynamics. Actually, due to the concentration of con- 
trol elements around hubs, the fact that a large number 



Marta Balbas Gambra, Erwin Frey: Social dynamics with peer support on heterogeneous networks 



9 




of small nodes are left unpoliced plays a more important 
role in the time evolution of the system than blocking the 
hubs. A large number of not so well connected agents are 
free to influence the society and drive it to a stationary 
state in which mafias do not get extinct. 

The micro-phases induced by the presence of control 
elements are now localized around highly and sparsely 
connected nodes depending on the police distribution. We 
have observed that policed areas are quickly occupied by 
citizens during the initial stages of the dynamics, while at 
later times the dynamics evolves in the unpoliced regions. 
In the stationary state, there is also a clear division in 
the degree of the nodes occupied by citizens and mafiosi: 
the population distribution follows the police distribution 
with citizens clustering around policed edges; see Fig. [TU] 
for an example with control elements drawn from the qi- 
distribution. 
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Figure 10. [color online] Stationary fraction of nodes of degree 
k populated by citizens (free squares) and mafiosi (red circles) 
for the case where control elements preferentially attach to 
hubs (qi). The simulations were run for a = 10, p = 0.2 in a 
7 = 2.5 scale-free network of size N = 8000. 



What is the system's behavior when control elements 
attach to specific targets as compared to the random case? 
The relation between the fraction of edges accommodat- 
ing control elements and the fraction of nodes which are 
protected by those edges seems to be the key factor in the 
analysis of the problem. Figures [TT] and [T^] show the frac- 
tion of both populations and the extinction probability for 
the three distributions, q±, q2, and random, as a function 
of network homogeneity, 7. 

The resulting fraction of controlled nodes, if control 
elements attach to large nodes, q±, whose degree is close 
to the maximal value Um , is smaller than in a random dis- 
tribution, p/k>M < p/ (k). Similarly, the effective fraction 
of protected nodes is larger than in the random case if 
control elements target sparsely connected nodes with the 
minimal degree k m , p/k m > pj (k). Thus, for a fixed police 
fraction, the population of mafiosi in equilibrium is larger 
for the gi-distribution than for the random case, and this 
population is larger for the latter than for the 92 distribu- 
tion: m(qi) > m(rand) > mfa), as illustrated in Fig. fTTl 
The difference in the stationary population fractions in- 



a 

o 




Figure 11. Stationary population fractions for different dis- 
tributions of the control elements and fixed police fraction 
p = 0.1 as a function of the increasing homogeneity of the 
scale- free network (7). The solid symbols represent the mafia 
fraction, the empty ones the citizen population in the station- 
ary state. The average degrees with increasing 7 ticks read 
(k) = 7.63, 4.31, 3.16, 2.69, 3.45. 



duced by the specific distributions naturally increases with 
the heterogeneity of networks (small 7) . 
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Figure 12. The police fraction needed to achieve extinction 
of mafiosi on a 7 = 2.5 ((/c) = 4.31) scale-free network of size 
N — 10000 crucially depends on the distribution of control 
elements on the network. A larger mafia fraction is needed 
when control elements are located around hubs (qi) compared 
with the random distribution, whereas the required fraction to 
fight the mafia depends on the initial conditions for control 
elements attached to small nodes. Initial conditions are mo = 
co = 1/2. 



Fig. IT2l illustrates that the fraction of control elements 
at which extinction takes place increases if the control 
elements concentrate around hubs (qi). In this case, an 
effective homogeneous structure is left unpoliced and, as 
we have discussed before, even small fractions of mafiosi 
manage to invade a citizens' population. In contrast, if the 
control elements preferentially surround small nodes ((72)) 
the unpoliced region consists of highly connected nodes, 
whose dynamics is similar to that of the well-mixed popu- 
lation for the unpoliced asymmetric model. In this case, 
the stationary state critically depends on the initial con- 
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ditions, so that we cannot draw any general conclusions 
on the effect of the distribution on the system's behavior. 
Whether a larger or smaller fraction of control elements 
than for the random distribution is needed to drive the 
mafia to extinction depends thus on the particular initial 
conditions in the unpoliced areas. 

6 Conclusions and outlook 

We have proposed a novel way to model the influence of 
the neighborhood in social dynamics which accounts for 
the protection of alike peers proportional to their frac- 
tion. In this way, both coexistence of opinions and con- 
sensus on the neutral position (citizens, in the language 
used in the introduction) emerge as possible solutions in 
a natural way. We find that the dynamics of the model 
in structured societies drastically differs from that pre- 
dicted for well-mixed populations on complete graphs. In 
the parameter window where the population dynamics ex- 
hibits bistability on complete graphs, we find coexistence 
between citizens and mafiosi on scale-free networks (and 
square lattices) . This is traced back to the local character 
of the interaction which allows mafiosi to invade societies 
in structured societies with small neighborhood sizes. We 
have rationalized this effect employing both heuristic ar- 
guments on the population dynamics as well as an exten- 
ded mean-field analysis, which accounts for the local net- 
work structure. Moreover, stochastic simulations for scale- 
free networks with increasing average degree (k) show that 
the coexistence region shrinks and the bistable regime is 
recovered in the limit of a complete graph |35| . 

The presence of control elements (police) splits up the 
social structure into unconnected subnetworks ruled by 
different dynamics. Neutral individuals (citizens) populate 
regions influenced by the control elements (policed areas), 
while the remaining (unpoliced) areas develop the same 
dynamics as in the model without police, yielding a sta- 
tionary state which depends on the network homogeneity. 
We have discussed the positive effect of this micro-phase 
separation on the population size of the minority group. 
In addition, we report a markedly sharp crossover from 
extinction of the minority to coexistence of species as a 
function of the heterogeneity of the underlying scale-free 
network. By preferential attachment of the control ele- 
ments to highly (sparsely) connected nodes, smaller (lar- 
ger) effective fractions of the control elements emerge. De- 
pending on the kind of nodes the control elements are pref- 
erentially attached to, the areas free of external influence 
(no police) exhibit now a degree of heterogeneity differ- 
ent from that of the whole network. We have specifically 
investigated the role of two different distributions of the 
external elements on the given structure, showing that 
the attachment rules for the control elements determine 
the threshold fraction needed to bring the minority group 
(mafia) to extinction. 

We conclude with a comparison of the mafia model 
with other non-linear opinion model and generalizations 
thereof. The mafia model can be viewed as an extension 
of the Abrams-Strogatz (AS) model [35] in the following 



sense. In the limit d — > and <fi — ¥ 0, where there are 
no more empty places in the mafia model, it reduces to 
a two-state reaction scheme c m with transtion rates 
w cm = s m m 2 and oj„ ic = s c c 2 , i.e. an AS model with a 
low volatility of a = 2. In the limit of the asymmetric 
model, s c — > 0, considered in this manuscript, the only 
stable fixed point is (m, c) = (1, 0), i.e. consensus on one of 
the two opinions as in the AS model. For general prestige 
values, s m / c ^ 0, the coexistence fixed point is unstable 
and both consensus states become stable, i.e. the dynamics 
exhibits bistability. 

These results highlight the importance of empty places 
(a third state) for the observed dynamics in the mafia 
model. Note, however, that the presence of empty places 
is only a necessary but not a sufficient condition for the 
emergence of bistability and coexistence. The functional 
form of the transition rates, i.e. the fact that citizens pro- 
tect themselves, also matters in an important and inter- 
esting way. The transition rates of our model, e.g. uj cm = 
s m m(l — c), weaken the tendency to switch to the opposite 
state by a factor (1 — c) as compared to neutral voter- like 
models. Instead of these transition rates, we could have 
adopted an AS scheme to investigate the effect of peer 
support: Lu crn = s m m a with 1 < a < 2. However, these 
non-linear transition rates are qualitatively different from 
the mafia model: for a fixed fraction of empty sites, 0, the 
rates in the mafia model are linear or non-linear depend- 
ing on the relative fractions of species and empty sites, 
while they are always non-linear in the AS model. 

Indeed, if one adopts transition rates like those of the 
AS model uj cm = cr m m a and uj mc = 0, i.e. generalizes 
the AS model to contain empty places, one does not find 
all the regimes of the mafia model. For a = 1, there is 
a transition between two monostable regimes, extinction 
m = and coexistence at a = 1/(1 + /3). For a = 2 
there is a saddle-node bifurcation with a transition from a 
monostable regime with extinction of mafiosi to a bistable 
regime, where extinction and coexistence are stable fixed 
points. However, it does not show a monostable regime 
with coexistence of both species. The mafia model exhibits 
richer behavior than its corresponding limits in terms of 
the generalized AS model, i.e. a = 1 and a = 2. In a 
sense, its dynamics combines all the regimes found for 
the two limiting cases of the AS model with empty sites: 
monostable regimes, with coexistence and extinction of 
mafiosi, and a bistable regime. 

Our model still lacks the inclusion of the evolution of 
the network of contacts which takes place in actual societ- 
ies. Considering adaptive networks which also account for 
topological features is, in our view, the most compelling 
extension of the model. 
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A Local mean-field approximation 

A standard mean-field description as given by Eqs. (|3J|- 
(7J) does not properly capture the phase diagram of struc- 
tured societes: the bistable regime of the mean-field the- 
ory largely becomes a coexistence regime in structured 
societies. For a well-mixed population, a standard mean- 
field theory which describes the interaction of every agent 
with all others is expected to correctly predict the sta- 
tionary state of the system. However, we find from nu- 
merical simulations that in structured societies (networks 
and square lattices) with a one-to-all interaction scheme 
this standard mean-field behavior is not recovered even 
in the limit where individuals diffuse very strongly [35] . 
The main reason is the one-to-all character of local in- 
teractions: while averaged frequencies properly reproduce 
the essence of pair interactions between two individuals 
(one-to-one), as in the voter model [TT) or in the rock- 
paper-scissors game [JT], it fails to capture the dynamics 
of one-to-all interactions, due to the discrete character of 
the neighborhood's composition. 

In the following we illustrate the relevance of the type 
of interaction for the phase behavior for square lattices 
without control elements. Here, the composition of the von 
Neumann neighborhood cannot take arbitrary values, but 
instead the possible fractions for both species and empty 
places are given by 0, 1/4, 2/4, 3/4, 4/4 with the con- 
straint that c + m + <fi = 1. There are thus fifteen possible 
neighborhoods Q t = (C 4 ,M 4 ,^) = (4,0,0), (3,1,0), 
(0, 0, 4), with Ci, Mi, and <Pi denoting the total number of 
individuals of each species (in a given neighborhood) ran- 
ging from to 4. Consequently, the transition rates U) mc 
and u cm are no longer continuous in structured societies. 

The number of possible compositions for a neighbor- 
hood is finite in networks and lattices. The smaller the 
neighborhood is, the smaller the number of possible com- 
positions and thus the further it is from the asymptotically 
continuous case of the global population. To account for 
the finite number of available neighborhoods one has to 
replace the averaged values of the species' frequencies en- 
tering the transition rates, Eqs. (JTJ) and ([2]), in a mean-field 
description by a sum over all possible neighborhoods f2i 
with their corresponding population fractions, weighted 
by the probability pi for every neighborhood^]: 



Won = crm 



(l-c)^5>cr^i(l-^] 



Wrnc — CTC 



(1 - m) — > J2 p * a T ( 1 
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Mi 
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(14) 



of drawing the specific combination out of a mixed sample 
with Nc citizens, Nm mafiosi and N<f> empty places. The 
corresponding probability for a given neighborhood is thus 
given by the multinomial distribution: 



Pi =p(d,Mi,$i) 



{Cj + Mj+^y. 
CAMA&A 



<j>*'. (15) 



Note that for a square lattice Ci+Mi+<Pi = 4. Accounting 
for these probability distributions of possible neighbor- 
hoods in the generalized transition rates, Eqs. (|13[) and 
(fl"4")l . we find that the fixed points of the ensuing exten- 
ded mean-field theory are still an absorbing and two co- 
existence states as for the standard mean-field, but the 
functional form of the coexistence solutions now read: 



0, 



(-2-/3-/3 2 )a± v /(l + /3)%( g -3) 
3a(l + /3) 2 



(16) 
(17) 



The lines at which the saddle-node, er* n = 3, and tran- 
scritical bifurcations, a tc = 4(1 + /3) 2 /(4/3 + /3 2 ), occur are 
drastically shifted with respect to those for the standard 
mean-field theory. The parameter region of the bistable 
regime in largely reduced to a thin stripe, and the coexist- 
ence phase now covers a much larger area in the parameter 
space, cf. Fig. [T31 
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Figure 13. Stability diagram for the asymmetric model in 
terms of the birth rate ft and the strength a within an ex- 
tended mean-field approach. Note that the bistable region has 
largely shrunk as compared to the standard mean fief shown 
in Fig. [2] The stability of the absorbing state m* = is lost 
and, accordingly, the parameter space for coexistence of both 
species is strongly enlarged. 



In our extended mean-field approach we still assume that 
there are no site correlations. Therefore, the probability to 
have Ci citizens, Mi mafiosi, and <I>i empty places is that 

4 To generalize the procedure for scale-free networks, one has 
to take the sum over the possibie neighborhood's composition, 
after having considered a sum over the possibie sizes of the 
vicinities, i.e. the degree k, weighted with the degree distribu- 
tion. 



The local mean-field analysis discussed here nicely cap- 
tures the dynamics for structured populations, accounting 
for the loss of stability of the absorbing state for a large 
region of the parameter space in favour of coexistence in 
structures. The one-to-all interactions turns out to be a 
key feature of the mafia model to understand its reach 
dynamics, and in particular the hegemony of species co- 
existence in a variety of structures. 
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